INTEGRAL TRANSFORMATION AND DARBOUX 
TRANSFORMATION OF HEUN'S DIFFERENTIAL EQUATION 
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Abstract. We review Darboux-Crum transformation of Heun's differential equa- 
tion. By rewriting an integral transformation of Heun's differential equation into a 
form of elliptic functions, we see that the integral representation is a generalization 
of Darboux-Crum transformation. We also consider conservation of monodromy 
with respect to the transformations. 



1. Introduction 
Heun's differential equation is defined by 

\ f d\ 2 f'j 5 e \ d af3z — q \ 

(1) \{Tz) + [- Z + — + —t)d- Z + z{z-l){z-t)) V = ^ 

with the condition 7 + 5 + e = a + j3 + l ( [TUl |TT] ) . It has four singularities {0, l,t, oo} 
on the Riemann sphere C U {oo} and all the singularities are regular. Heun's dif- 
ferential equation is a standard form of a second-order Fuchsian differential equation 
with four singularities. Heun's differential equation and the differential equations of 
its confluent type appear in several systems of physics including analysis of black 
holes, AdS/CFT correspondence ([7]), crystal transition ([II]). fluid dynamics ([2])- 
A standard form of a second-order Fuchsian differential equation with three singular- 
ities is the hypergeometric differential equation, which is celebrated in physics and 
mathematics. Although analysis of Heun's differential equation is much difficult than 
that of the hypergeometric equation, several solutions are known. For example, Heun 
polynomials ([10]) are polynomial solutions of Heun's differential equation, which are 
related with quasi-exact solvability. Finite-gap integration is applicable to Heun's 
differential equation for the case 7, 5, e, a — [3 G Z + l/2 ([231 HI 021 EES]), and mon- 
odromy of Heun's differential equation can be calculated in terms of hyperelliptic 
integral Q15J) and by Hermite-Kricherver Ansats ([16]). 

We now change variables of Heun's differential equation. Let p(x) be the Weier- 
strass doubly-elliptic function, which has double-periodicity p(x) = p(x + 2ui) = 
p(x + 2u 3 ). We set ujq = 0, uj 2 = —lji — lu 3 , = p(uii) (i = 1, 2, 3) and 

(2) z = eg^, t=^, y*(z) = f(x), *{z)=z=} L (z-l)=£(z-t)=£. 
Then Heun's differential equation ([1]) is transformed to 

(3) H (io,h,hM f{x) = H <h,hM>) = + k(k + l)p(x + Ui), 
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where 



(4) l = P-a- 1/2, h = - 7 + 1/2, Z 2 = -5 + 1/2, h = -e + 1/2. 

The singularities {0,1, t, oo} of Eq.flT]) correspond to half-periods {u 1 ,u 2 ,uj 3 ,u } of 
the elliptic function p(x). The parameter q essentially corresponds to the eigenvalue 
E (see [H]). Note that the operator JjvoMiM) is Hamiltonian of Inozemtsev system 
of type BC\. Lame's differential equation is obtained as a special case l\ = 1% — I3 = 
( 1 = 5 = e = l/2). 

By the way, Darboux transformation is extensively applied to soliton theory ([9]), 
and one of the features is isospectrality. Here we briefly review Darboux transforma- 
tion for Schrodinger equations, which is simpler than that for soliton equations. 

Set H = — + q(x) and assume that 4>o(x) is an eigenfunction, i.e. H(f)o(x) = 
Eo<fio(x) for some E . Then q(x) = c/)q(x) /4>o{x) +E . Define an annihilation operator 
and a creation operator by L = £ - gjg, V = -£ - gg. Then L*L = -£ + 

+ (Si) 2 = H - E * ^ ^changing and L, we have W = -£ - 

(jfoffi) + (|om) = H — 2 ^gfj J — £?o- Hence we obtain the Schrodinger operator 

H = + + 2 (jgg))' and the relations H = L^L + Eq, H = W + E . It 

follows from the relation LH = LL^L + E L = HL that, if f(x) is an eigenfunction 
of H with the eigenvalue then Lf(x) is an eigenfunction of H with the eigenvalue 
E, because H(Lf(x)) = LHf(x) = L(Ef(x)) = E(Lf(x)). The transformation from 
H to H (or the operation by L) ia called Darboux transformation. Note that the 
operator L annihilates the 1-dimensional space C<f)o(x). 

Darboux transformation is applied to isomonodromic property of Heun's differen- 
tial equations, which we will explain by an example. Let if ( 2 > '°>°) be the operator 
defined in Eq.((3]). It follows from a straightforward calculation that the function 
4>o(x) = (p(x) — e\)(p(x) — e-i) is an eigenfunction of _£/"( 2 ' > >°) with the eigenvalue 
3e3, i.e. H( 2 ' 0,0,0 >(j)o(x) = 3es4>o(x). The annihilation operator and the creation oper- 
ator arp PYnrpwH fl c T — A — p'W p'W rt — A _ 

ator are expressed as L - dx 2(p(s) _ ei) 2 ( P (x)-ea)' ^ ~ <fc 2(„( ai )-e 1 ) 2( P (*)-ea) 

and we have H = H + 2 (^|y)' = fft 1 ' 1 ' 1 ' ), H^^L = L#( 2 .oao) The operator L 

determines a map from an eigenfunction of _£/~( 2 >°> >°) with an eigenvalue E to the eigen- 
function of H ( 2 '°'°'°) with the eigenvalue E for all i?. Since the operator L is doubly- 
periodic, the monodromy of solutions to (_£/~( 2 >°>o,o) —E)f(z) =0 coincides with the one 

to (H^)-E)f(z) = 0. In other words, the operator L = 

induces isomonodromic property of _£/"( 2 >°> ' ) and H^ 1 ' 1 ' 1 ' '. Note that monodromy 
of solutions to (#( 2 <°>o,o) _ E)f(z) = (resp. (H^ 1 ^ - E)f(z) = 0) had been 
calculated individually by finite-gap integration ([16J), and it was found later that 
Darboux transformation provides a direct connection of coincidence of monodromy 
([HIELTI)- Calculation of monodromy is important for application to physics, because 
monodromy is used judge periodicity of eigenf unctions and periodicity is related to 
boundary condition for the eigenf unctions. In this paper we generalize Darboux trans- 
formation and apply it to isomonodromic property for Heun's differential equation 
by following |17j . 
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By the way, an integral transformation on Heun's differential equation was noticed 
in a connection with middle convolution ([19], EI]), although it had been established 
before by Kazakov and Slavyanov [5] . In this paper we will study the integral trans- 
formation on Heun's differential equation and observe relationship with generalized 
Darboux transformation. 



2. DARBOUX-CRUM TRANSFORMATION 

We introduce a proposition that is a generalization of Darboux transformation. 

Proposition 1. M, W Let H = + q(x) and U be a n-dimensional space of 
functions which is invariant under the action of H. Let 

(5) l=(i)* + Sa«(')(ir 

be the operator that annihilates any elements in U , i.e. Lf(x) = for all f(x) G U . 
Set H = —-^2 + q{x) + 2c[(x). then we have 

(6) EL = LH. 

We call L the generalized Darboux transformation or Darboux-Crum transforma- 
tion. Note that Crum [3] had obtained a proposition whose expression is slightly 
different. If n = 1, then we reproduce Darboux transformation. 

In order to apply Proposition [T] for Heun's differential equation, we recall quasi- 
solvability of Heun's equation. 

Proposition 2. ( |14[ Proposition 5.1}) Let ojj be a number such that oti = — or 
Oii — h + 1 for each i G {0, 1, 2, 3}. Set d = — Y^=o and assume d G Z> . Let 
K*o,au,a2,a3 be the d + 1-dimensional space spanned by I $(p(x))p(x) n } , where 

' ' ' (, J n=0,. ..,d 

= (z — ei) ai / 2 (z — e 2 ) a2 ^ 2 (z — e 3 ) a3 / 2 . Then the operator JJ"(Wi>W3) preserves 
the space V^, Q , 1jQ , 2jQ , 3 . 

Let £ Q0jai)Q2)Q3 be the monic differential operator of order d + 1 which annihilates 
the space V a0 ^ aita2)Olz . By calculating the coefficient c\(x) in Eq.([5]) for the operator 
-^«o,«i,«2,Q3 5 we have Darboux-Crum transformation for Heun's differential equation. 

Theorem 3. fl!7[ Theorem 3.3],) Let c*i be a number such that oti = —li or ai = l{ + 1 
for each % G {0, 1, 2, 3}. Set d = — Y2^=o a i/^ an d assume d G Z> . Then we have 

(J\ fr(a +d,ot-i+d,a 2 +d,a 3 +d) r T fj(lo,h,h,h) 

Let fi(x,E), f 2 (x,E) be a basis of solutions to (ifCWi.Ws) -E)f(x) = 0. Since the 
operator ff(Wi,W3) is doubly-periodic, the functions + 2ujk,E), f 2 (x + 2uk,E) 
(k = 1,3) are also solutions to the differential equation. Let M 2LUk (E) (k = 1,3) be 
a monodromy matrix on the shift x — > x + 2^ with respect to the basis £"), 
/ 2 (x,£)}, i.e. 

(8) (f 1 (x + 2u k7 E) f 2 (x + 2uj k ,E)) = (ft(x,E) f 2 (x, E))M 2u>k (E). 

It follows from absence of first order derivative of the differential equation that 
det M 2 ^ fe (-E) = 1 {k = 1,3). Note that tiM 2u , k {E) is independent from the choices of 
a basis of solutions to (iyPo.'i.'a.fe) _ E)f(x) = 0. Set fi(x,E) = L ao ^ ua ^ a Ji{x,E) 
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(i = 1,2). Then H( ao+d > ai+d > aa+d > aa+d )f i (x,E) = Ef t (x,E). Since L 
doubly-periodic, we have 

(9) (A(x + 2u k , E) } 2 {x + 2u k , E)) = E) f 2 (x, E))M 2uJk (E). 

Hence the monodromy structure of H^ l °' 11 ' 12 ' 1 '^ with respect to a shift of a period 
coincides with the one of H ( - ao+d ' ai+d ' a2+d,a3+d \ and the operator L ao aijOC2jOC3 defines 
an isomonodromic transformation from HfoMMa) to jj{a +d,ai+d,a 2 +d,a 3 +d) _ j n p ar _ 
ticular, we have tiM 2uJk (E) = trM 2uJk (E), where trM 2uJk (E) (resp. trM 2uJk (E)) is 
the trace of the monodromy matrix of solutions to (_ff"(Wi,W3] _ E)f(x) = (resp. 

fjj((x +d,ai+d,a 2 +d,a 3 +d) _ — 0) 

We now provide an example of Darboux-Crum transformation. On the case Iq — 2/ 
(/ G Z>i), h = h = h = 0, we set «o = —2/, «i = a 2 = 1, «3 = 0. Then 
d = — (cko + • ■ ■ + ^3) /2 = / — 1 and we may adopt Theorem [31 We have 

(10) iftW'^L^^Q = L_ 2Z , M , tf^ A0) ; 

because H^~ l ~ l ' l,l ' l ~ l > = H^ l,l,l ' l ~ l \ Hence _£/"( 2Z > > >°) is isomonodromic to ijW'M-i) f or 
/ G Z> x . If I = 1, then d = 0, H^ l ^L_ 2 ^ lfl = L_ 2ilili0 # (2 ' '°' 0) , and the operator 
£-2,1,1,0 is written as £-2,1,1,0 = £ - 2 J^ ei) - 2(p $l 2 ) - Namely we reproduce the 
example in introduction. 

3. Application to finite-gap integration 

The method of finite-gap integration has been studied more than 30 years (p3|) and 
a target of earlier study was KdV equation. It was also applied to Heun's differential 
equation ([231 IU E2l EEE]) with the condition / , Z 1; l 2 , Z 3 G Z and to a certain class of 
Fuchsian differential equations ([20]). Consequently we have formulas of monodnomy 
of Heun's differential equation by hyperelliptic integrals ([15]), Bethe Ansats Q13J) 
and Hermite-Krichever Ansatz ([16]). A remarkable feature of finite-gap integration 
is existence of an odd-order differential operator A which commutes with the operator 

Jj{h,h,hh) ^ [ e ft (Wi>W3)l — Q 

We can construct the odd-order differential operator A by composing Darboux- 
Crum transformations for the case Iq, h, h, h G Z. 

Proposition 4. f|17j ) If lo, l\, l 2 , I3 G Z, then we can construct an odd-order dif- 
ferential operator A such that [A, /f(Wi>W3)] = by composing four Darboux-Crum 
transformations. 

For details, see [TTJ, Proposition 6.1]. If lo = 2, li = l 2 = l 3 = 0, then there exists a 
differential operator A of fifth order written as A = £2,-1,-1,0-^1, -2, 1,0-^0,2, -1,-1-^-2,0,0,0 
which satisfies [A, H^ ' ^} = 0. 

4. Integral transformation of Heun's differential equation 

We introduce integral representation of Heun's differential equation as a general- 
ization of Darboux transformation. 

Integral transformation on Heun's differential equation has been noticed in a con- 
nection with middle convolution ([121 EI]), although it had been established before 
by Kazakov and Slavyanov [5]. In order to formulate the integral transformation, we 
fix a base point o of the integrals in the complex plane C appropriately. Let p be an 
element of the Riemann sphere C U {00} and 7 P be a cycle in the Riemann sphere 
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with the variable w which starts from w = o, turns the point w = p anti-clockwise 
and ends at w — o. Let [j z , j p ] = lzl P lz l % 1 be the Pochhammer contour. 

Proposition 5. fl5l EE]j (i) Set 

{// - (2 - «)}{// - (2 - /?)} = 0, i = 7 + fi - 1, 8' = 5 + n - 1, e' = e + fi - 1, 

(11) a' = //, = 2/x + a + /? - 3, q' = q + (1 - /j)(e + St + (7 - + !))• 
Lei y(w) 6e a solution to 

(12) S + fH^ + ^r + TnirTi^O' 

> « aao^ \m «)— 1 w—tj aw w(w—l)(w—t) t ' 

Then the functions (i G {0, 1, t, 00}) 

(!3) = I bzm] y(w)(z - w)-»dw 

are solutions to 

(14) ^ 4- 4- 4- ^ 4- ?; - 

(ii) Under the notation of (i), we assume fi G Z>i additionally. Then 

(is) yW = (i) / " 1 y(«) 

is a solution to Eg. HM. 

We rewrite Proposition [5] to the form of the elliptical representation by setting 
fi — d 4- 2. It is remarkable that the eigenvalue E is unchanged by the integral 
transformation. 

Proposition 6. f|22] ) (i) Let <r{x) be the Weierstrass sigma function, <Ti(x) (i = 
1, 2, 3) be the Weierstrass co-sigma function which has a zero at x = Ui, and Ii (i = 
0, 1, 2, 3) be the cycle on the complex plane with the variable y such that points y = x 
andy = —x+2uji are contained and the half-periods Zui+Zus are not contained inside 
the cycle. Let a% be a number such that a% = —k or ct{ = U + 1 for each i G {0,1,2,3}. 
Set d = -ELo^/ 2 - Let f( x ) be a solution to (#0°>'i>W 3 ) _ E)f(x) = 0. Then the 
functions 

(16) f(x) = a(x) ao+d+1 a 1 (x) ai+d+1 a 2 (x) a2+d+1 a 3 (x) a3+d+1 ■ 

I h f(y)a(y) 1 - ao a 1 (y) 1 -^a 2 (y) 1 -^a 3 (y) 1 ^(a(x + y)a(x - y))- d ' 2 dy 

(i G {0, 1,2,3}) are solutions to (H( ao+d > ai+d > aa+d > aa+d > _ E)f(x) = 0. 

(ii) Under the notation of (i), we assume d G Z>_i additionally. Let f(x) be a 
solution to (ifCo'^'^'^) _ E)f(x) = 0. Then the function 



(17) /» = p'{x) d+1 f[(p(x) - e^ 2 (^y^) ^ \m {[(Pi*) - e 

is a solution to (#(ao+d,ai+d,a 2 +<f,a 3 +d) - E)f(x) = 0. 

If d G Z> and f(x) G V a0taita2tO(S , then the function f(x) in Eq. ljTTI) is identically 
equal to zero, and it follows that 



( oti/2 



P'(x) d+1 



3 / 1 d \ d+1 3 

Y[(p(x) - ei) a ' /2 o —T-r-r- o Y[(p(x) - e l )~^ /2 = /-,„,,„.», 
i=i \P{X)ax/ , =1 
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Therefore we may regard the integral transformation as a generalization of Darboux- 
Crum transformation by removing the condition d G Z> . 

We can express the monodromy of solutions to i^}j ( ' a(>+d,ai+d ' a2+d ' a3+d ^ — E)f(x) = 
in terms of the monodromy of solutions to ^H^ ' ll > l2 ' 13 ' — E)f(x) = by applying 
integral transformations, which was performed in [22J. 

Theorem 7. ( |22j ) Let at be a number such that oti = —l{ or ati = k + 1 for each i G 
{0,1,2,3}. Set d = -J2 3 i= o a i/ 2 - Let he {1,3} and M 2uJk (E) (resp. M 2u>k (E)) be the 
monodromy matrix by the shift of the period x — > x+2uk with respect to a certain basis 
of solutions to (ifCWi.Ws) _ E)f(x) = (resp. (#to+«*.«i+d,«a-H*,«3+«0 - E) f(x) = 0). 

Then tiM 2uJk (E) = tiM 2ulk (E). 

Hence the monodromy structure of ij(Wi,W3) w ith respect to a shift of a period co- 
incides with the one of H^ 0+d ' ai+d ' a2+d ' a3+d ^> for all d(= - Yf i=0 an/2). As a corollary 
we have 

Corollary 8. ([22\) We keep the notations in Theorem^ Let k G {1,3}. If 
there exists a non-zero solution f(x,E) to (_H"(Wi>W3) _ E)f(x,E) = such that 
f(x + 2uk,E) = Ck{E)f(x,E), then there exists a non-zero solution f(x,E) to 
( H (a +d,a 1+ d,a2+d,a 3 +d) _ E)f(x,E) = such that f(x + 2u k , E) = C k (E)f(x,E). 
In other word, periodicity is preserved by the integral transformation. 

For example the monodromy structure of H^ 21 ' ' ' ^ with respect to a shift of a 
period coincides with the one of H^ 1,1,1,1 " 1 ^ for all I G R. In particular, the condition 
of eigenvalues E such that there exists a non-zero periodic solution to H^ 21 ' ' ' ^ f(x) = 
Ef(x) with respect to the shift x x + 2uj\ coincides with that of eigenvalues E such 
that there exists a non-zero periodic solution to H^ 1,1,1,1 " 1 ^ f{x) = Ef(x) with respect 
to the shift x — > x + 2u>i for all / G R. 

5. Concluding remarks 

In this paper we observed some aspects of Darboux-Crum transformation and in- 
tegral transformation on Heun's differential equation. In particular we can regard 
Darboux-Crum transformation as a specialization of integral transformation, and 
periodicities of solutions to Heun's differential equation is preserved by the transfor- 
mations. We may interpret our result as a different realization of the observation by 
Khare and Sukhatme [6]. 

We note that integral transformation may produce new solutions to Heun's differ- 
ential equation. For example, solutions of Heun's differential equation with the con- 
dition 7, 5, e, a + 1/2, (3 + 1/2 G Z (i.e. l , h, h, h^Z+ 1/2, l + h + 1 2 + 1 3 G 2Z + 1) 
can be expressed in use of finite-gap solutions, and the special case 7 = 5 = e = 1, 
a = 1/2, (3 = 3/2 (i.e. l = 1/2, l x = l 2 = l 3 = —1/2) was previously studied by 
Valent [23] to understand an eigenvalue problem related to certain birth and death 
processes. We hope further applications of new solutions to physics and mathematics. 
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